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I. INTRODUCTION 



Symmetries of partial differential equations are used for description of the general set of 
solutions, for producing families of solutions from known exact solutions, for construction 
of invariant solutions etc-*^"^. 

Intcgrablc evolution equations (systems of equation), with the group point of view, has a 
remarkable property: they come in hierarchies of commuting flows which often it is possible 
to describe with the help of suitable recursion operator as 

= A-^Ki, 

where A is recursion operator and Ki is some vector field. In many cases Ki[u] = u^^, 
which presents a shift of spatial variable x G R^. The fact of compatibility of these flows 
make possible to say about simultaneous solution u{x, t2, ts, ...) and investigate their group 
properties. 

In this article we investigate integrable evolution equations (systems of equations) asso- 
ciated with stationary Schrodinger equation with energy-dependent potential. In this case, 
recursion operator A is known in its explicit form. That allows us to establish transforma- 
tion property of A. We construct one-parameter linear point transformation generated by 
the shift of spectral parameter A — A = A — e. Particular case of these transformations 
is well known Galilei transformation for Korteweg — de Vries equation. Also we construct 
scaling invariance group. We next apply standard symmetry method to construct simulta- 
neous solutions, invariant w.r.t. Galilei-like group QN,n- Corresponding ansatz involve the 
vector-function U = U(X, T2, ...,Tn_2,Tn), where {X = Ti,T2, ...,Tn-2,Tn} is the set of 
differential invariants. We show that this ansatz reduce first N members of hierarchy 

Ut^=A"*-i[u]u, = K^[u], m = l,...,A^ 

to (A^ — 2) members of that, i.e., 

Ut^ = K^[U], m = l,...,A^-2, 

with additional Galilean and scaling self-similarity constraint. 

The plan of this article is as follows. In Sec. II, we construct symmetry transformation 
groups. Sec. Ill is devoted to construction ^iv,n- invariant solutions. 



II. INTEGRABLE HIERARHIES AND SYMMETRY 
TRANSFORMATIONS 

Let us recall some relevant notions which are useful throughout this paper. Let M be 
a manifold of smooth vector-functions u : ^ C". We denote by A^, the algebra of 
polynomials in finite collection of variables u^^^ 

Definition 1: Let K[u] = (A'i[u], A'„[u])"^ G TuM is a vector field and A : TuM 
T^M is a linear operator. The Gateaux derivatives of K and A with respect to u in the 
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direction X G T^M are defined through the relations 



K'[u](X) = ^ 



T=0 



K[u + rX], A'[u](X)= ^ 



A[u + rXl. 



r=0 



as 



Definition 2: The Lie derivatives of K and A in the direction X are defined, respectively, 
LxK = [X, K] = K'(X) - X'(K), LxA = A'(X) - [K', A]. 



The linear space T^M endowed with the commutator [X, Y] = LxY bears the structure 
of the infinite Lie algebra Vect. 

Definition 3: The operator A : T^M T^M is called hereditary if its Nijenhuis torsion 
vanishes i.e. 

rA(X,Y) = [AX,AY]-A{[AX,Y] + [X,AY]} + A2[X,Y] = 
for any X, Y G T^M. 

Vector fields A'^~^Ki span abelian subalgebra in Vect if an operator A is hereditary and 
Lk,A = ^ 



It is known that, with each auxiliary linear equation 



i=l 



one can associate the hierarchies of commuting fiows^ 

ut^ = K„[u] = A"~^[u]u^, u = 
with the hereditary recursion operator of the form 

/ ... + ui + ^ui^d-^ \ 



(1) 



(2) 



Aiul 



-1 
-1 



V 



M2 + |m2x<9^^ 

■•. : 

-1 Un + \Unxd^^ ) 



(3) 



It is obvious that L-^^ K = 0. 



In what follows, we define the operator occurring in (H) so that d~^{f) G A° for 
any f = f[u\ E Im CI A^, where C A^^ is a ring of differential polynomials in the 
fields Ui with zero constant terms. 
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A. Galilei-like symmetry transformations 

In the previous paper^ we have proposed a construction of Galilei-like symmetry trans- 
formations for equations (0). To make this paper self-contained, we gives all proofs below. 

Let us define the one-parameter linear transformation of the fields Ui by the shift of the 
spectral parameter A through the relation 



J2u,ix)i-Xy~' + (-A)" = Y^u,{x)i-X + e)-i + (-A + e) 

1=1 i=l 

Comparing the coefficients of the same powers of A, one obtains 

n 

ui = Fi(u, e) = 5]u,e'-^ + e", 



Ui+i = Fj+i(u, e) 



i=l 

1 d'Fi{u, e) 



(4) 



i = 1, n — 1. 



It is obvious that inverse transformation to (|^) is given by 

Mi=Fi(u, -e), Ui+i = Fj+i(u, -e), i = l,...,n-l. (5) 

Eq. @ can be rewritten in vector form 

u = F(u,e) =A(e)u + d(e), (6) 

where A{e) is an upper diagonal matrix with units in main diagonal and d(e) is a vector. 
From (3) we easy obtain 



Here and in what follows the symbol denotes binomial coefficient (^^^ 
Lemma 1: The recursion operator satisfy the following identity: 

A[u]|„=^(^,,)=A(e) iA[u] + e) A'\e) 

Proof: Let us denote 

pie) = po + ... + Pn-ie""' + e" = (^^^ + m + ^Ui,d~' 
By virtue of Eq. @ we can write 

/ ... 

-1 ... 



(7) 



(8) 



He) = A[u] 



U=F(U,€) 



ui=Fi{u,e) 



P{e) \ 



-1 



\ 



■•• 

-1 



1 r)("-l)| 

{n-l)\P ' 



e) ) 



It is obvious that A(0) = A[u]. Since the matrix A{e) is nondegenerate, we can rewrite Eq. 
(I) as 

A(e)(A(0) + e)=A(e)A(e). (9) 

The coefficients pi are operators but it is obvious that when proving fuffiUment of relation 
(P) we can operate with pi like with a numbers. 



We have 



A{e)A{0) 



— e 


-e2 .. 


gTl— 1 


-1 


-Cle .. 







-1 










-1 



p{e) — 



p'{e) 



eAie) 



( 



A(e)A(e) 



(e 6^ 
e 




V 




-1 -C\e 



\ 





^n— 1 


e" 




... o„_2e 


fil n-1 


e 


(^2 .ra-3 
... 0„_2t 










e 
















/^l ,n-3 



V 











(n-l) 



,p"{e) - cr 



(n-l), 



w-2. 



Comparing diagonal-wise matrix elements and using the relation 

^k-1 — ^k-1 5 

one can verify that (|^) is fulfilled. □ 

The relation (|) defines transformation law of recursion operator A[u] w.r.t. but 
it should be noted that operator d~^, occurring in A[u], now acts on a ring of differen- 
tial polynomials in the transformed fields Ui and its action must be consistent with this 
transformation. Let F,,, : Au — > and F^^ : A^ — >■ are maps of rings of differential 
polynomials generating by Then we must to require that 



1-1 



(10) 
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for any / [u] G Im dx C A^. 

So, when computing vector fields A'~[u]u2,. we must to take into account the condition 
(p^). For example, 

dx^i^nx) =Un + ne (11) 
because F^^ (n„ + ne) = Un E ^4° . Taking into account (|lT]) we have 

A[u]ux = K2[u] + -eux, (12) 

To compute A^[u]u2; there is a need to calculate the nth component of the vector field K2[u]. 
We have K2[u] = \uxxx + ^uux £ Inic?^ in the case n = 1 and i^'2,n[u] = —Un-i,x + ^UnUnx G 
Im dx for k >2. Taking into account (p!OD, we obtain 

dx'mu]) = \uxx + lu'-le', (13) 

d-' (ir2,n[u]) = -Un^l + ^< + (^^(n - 1) - ^n^) (14) 
Using ([T|), (|1|) and (|14D we obtain: 

ft 

A^[n]ux = A[u]ir2[u] + -eA[u]u^ 

= K3[u] + Qn(r2 - 1) - ^n^^ + (k2[u] + ^eu^. 

n /I 1 \ 

= K3[n] + -eK2[u] + i^-n{n - 1) - -n^j e^u,. (15) 

Using lemma 1 and Eqs. (^^), (p^ ) we obtain evolution equations on transformed functions 

Ui{x,t2,t3) 

= (A[u] + e) u^. 
= K2[n]+ (^^ + l)eu,, (16) 



= (A[u] + e) 
= A'^[u]ux + 2eA[u]u^ + e^u^ 

K3[n] + ^eK2[u] + (in(n - 1) - ) e% 



+2e (^K2[u] + -euxj + e'ux 
K3[u] + ('S + 2) eK2[n] + ^ (i + 1) (i + 2 ) e^u.. (17) 
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Let us define two one-parameter linear transformations of independent variables 



X 



Boie) 



X 



1 (f + l)e 



X 
t2 



and 



/ 


X \ 




1 X \ 




h 






V 


h) 




\h 1 





V 



1 





(f + 2) 
1 









(j) 


) 





(19) 



It is easy to check that matrices B2{e) and -63(6) satisfy group relation B{e)B{ei) = B{e+ei). 
By straightforward calculations, it can be verified that: 

(i) system of evolution equation Up^ = K2[u] = A[u]u^, is invariant under transformations 

(m and (in, 



(ii) a pair of systems u^^ = K2[u] = A[u]uj. and Ujg = Kafu] = A^[u]ua, is invariant 
under transformations (^ and (pj|). 

We can undertake to generalize transformations ( [TB| ) and (|T^). Let us denote t = 
(x, ^2, ^3, ^Af)^- We suppose that matrix B^i^e) defining transformation of independent 
variables 

t = 5^(e)t (20) 
is upper diagonal matrix with units in main diagonal and is written in the form 

/ 1 62,16 ... 6iv,ie^-' \ 

1 ■•. ; 



BN{e) 



bN,N-l^ 

1 



Here bk^i are some coefficients to be determined. Necessary condition in order that (^ 
complete transformation (^ to define one-parameter group is the relation 

Siv(e)Siv(ei) = 5jv(e + ei), 
which element-wise can be written as 



r-l 



bm-\-r,m (^ ~l~ ^l) bm+r,in 



i=i J 

h f:^ h h ■ <rV''~* -1-7) c'^ 



r-l 



i=l 



Thus the relation 



Cth 



'r+k,k "r+k,k+i"k+i,k 



(21) 



(22) 
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must be fulfilled for each i = l,...,r — 1. The solution of Eq. (p^ is uniquely given by 

1 

^- 1=1 

Thus to derive coefficients bmi there is a need to calculate &m,m-i- In order to do that we 
observe that 

ut^ = Km[u\ +b 

171,171—1 

eK„_i[u] + O(e^). (23) 

Taking into account (|T^) we have 

= A"^-^[u]u, + (m - l)eA™-2[u]u, + 0(6^) 

= A'"-2[u] (k2[u] + |6u,.) + (m - l)6A™-2[u]u,. + 0{e') 

= A'"-2[u]K2[u] + Q + m - l) eA'"-2[u]u, + ©(e^) 



= K„[il] + Q + m - Ij eK„_i[u] + ©(e^). 
Here we also consider that A[u]K2[u] = Kafu] + O(e^). From (|^) it follows that 

1 ^ /n 



' • i=i ^ ^ 



(24) 



In what follows, we will denote the groups of one-parameter transformations (^, ( pOD 
by QN,n- The infinitesimal generators (vector fields) of QN,n can be calculated immediately 
to yield 



N-l 



where Un+i = 1- For convenience, here we denote ti = x. As can be checked, respective 
characteristics are given by 

where A'{0) and d'(0) denote a valuation of derivatives with respect to group parameter 
calculated at e = 0, being A'{0) = 1]"=/ ii^j+i^j and 



d'(0) 



/ \ 


V n / 



Here E^j denotes matrix with unit in place and zeros elsewhere. 
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Repeatedly applying recursion operator A[u] to Qn,^ we obtain infinite sequence of 
characteristics QTv.n- particular, once applying A[u] to Qn^u yields the characteristic 



N 



where J = diag(ra, n — 1, 1). This characteristic corresponds to scaling invariance group. 
Below we construct scaling invariance transformations making use a transformation prop- 
erty of recursion operator A[u]. 

B. Scaling symmetry transformations 

Let us investigate invariance properties of the systems (|^) w.r.t. scaling transformations. 



We denote pi = \dl + Ui + lui^d^ ^, Pi = Ui + \ui^d^ \ 



relation 



/ 


V ... 



... \ 

■•. 

/?V 



/ 

-1 

-1 



= 2,...,n. It is obvious that 
pi\ 

P2 



P3 







/ 

-1 

-1 

V ... 



... -I Pn J 

/ Z?'^ 







V 



... \ 

■•. 

pj 



(25) 



V 

... /?>i \ 
... /3'^-V2 

■•• ■ r-'P3 

■•. ; 

-1 PPn J 

where (3 is arbitrary number, is valid. In what follows, (3 will be group parameter. 
Let us introduce scaling transformation of variables x and uf 

Consequently Eq. (^) we comes to the relation 

AMIu=A(/.)n = MP) (/5A[u]) A-\(3), (26) 
where A{P) = diag(/5", ...,/3). Using the Eq. ([26|), we obtain equation 



^(n/2)+m-l^m-l^-|- 
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from which it follows that by scaling evolution parameter tm as tm = we arrive 

at 

Thus we have the following statement: 
Corollary 1: The system of equations 

with the recursion operator is invariant under scaling transformation 

x = /5"/\, = Ui = p'-^-W (27) 



One can assign, to each variable, scaling dimension [■], so that [x] = n/2, [tm] = {n/2) + 
m — 1 and [ui] = i — n — 1. 

III. INVARIANT SOLUTIONS 

Given Galilei-like symmetry group QN,n, corresponding one-parameter family of solu- 
tions of (I) is given by 

u(e,x,t2, ■■■,tN) = A{e)u{x,t2, ...,tAr) + d(e) (28) 

Next we are going, using standard symmetry method, to construct ^jv^„-invariant solutions 
of equations (0). To do this, we need to find the set of generators of ^^^.n-invariant functions 
algebra. 

Proposition 1: For the group Qn^u, there exist functionally independent differential 
invariants in the form: 

T,=t,+ E (-1)W^^4^, j = l,...,N-2, 

k=l ''N 



Tn — tN, 

"-i (29) 
= «i + II(-l)''/?i,fcWi+fcWn, j = 1> -^n- 1, 

k=l 

_ 2ntN-i 
^n-Wn+ ^2iV + r2-2)t^' 

where the coefficients aj^k one? jSj^k ore positive rational numbers to be determined. 

Proof. It is evident that T/v is an invariant. Consider Tj, for j = 1, — 3. We have 

VjV,n(Tj) = +j) 
, ■ , i\ tj+2tN-l . . I ^-.N-j-l ("^ , Ar iA tfyft^Ji ^ 
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tj+1 



t]\f 

+«j,2 ;2 + - + (-1) «j,Jv-j-i jv=FT 

''AT % 



Require that VAr,„(Tj) = 0. Collecting similar terms and equating it to zero, we obtain 
recurrence relations: 

tj+i ■■ (l+j) = (| + A^-l)a,-,i, 



(f + ^ - 2) a,-,iv-,-2 = (TV - j) + ^ - 1) «.-,iv-,-i. 
Separately, consider the case j = N — 2. Assuming that ^^N,n{TN-2) = 0, we obtain relation 

Viv,n(T7v-2) = (I + - 2) tN-i - 2q;jv-2,i (I + - 1) tjv-i = 



1^ ^ (f + ^-2)«, 



from which it follows 



So, all coefficients aj^k are uniqely defined by recurrence relations above. ^From these 
relations, we easy obtain 

1 /n \ 

<^j,k-— Tin o + ' + ^'~^ ' i = i>->^-i-2, 

A;!(f + 7V-l) i=i^^ ^ 

(N-j-l) ""Tj' (n . . , 

[N - i){N - 1 -l)\i^ + N -l] i=i ^2 



[N-j){N-3-l)\{l + N 

By simple calculations, it is easy to prove that Vn is an invariant of QN,n- Consider Vj, 
for j' = 1, n — 1. Requiring that 'yN,n{vj) = we obtain 

V7V,n(^^j) = -n + 2(3j^2Uj+2Un + ••• + (-l)""-'"^(n - j - n_j_iM„_i<"^~^ 

+ - i)/5j,n-j<"^} - + (j + i-Uj+sWn - (j + 2)(3j^2Uj+zul + ... 

Collecting similar terms we obtain recurrence relations: 

Uj+i : n/3j- 1 = J, 
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u 



+kUn ■ knPj^k = {j + k- l)Pj,k-i, k = 2, ...,n - j - 1, 



from which it follows that 

Pj,k = '^^^l'^ , k=l,...,n-j -1 

fj. n-1 
P3,n-] ^_ j j^l ^n-j ■ 

It is obvious that all invariants determined above are functionally independent and generate 
the algebra of ^^^.n-invariant functions. □ 

It should be noted that [Tj] = [tj] and [vj] = [uj]. Denote T = (Ti = X, T2, ...,T/v_2)"^ 
and t = {x,t2, tAr_2)"^. In what follows, it will be useful following two lemmas. 

Lemma 2: The relation 

T = 5A._2(r)t + c(r), (30) 

where 

tN-l 

T - 



(f + iV-l)t 

and c{t) = (ci(r), CAr_2(T))'^, where 



N 



CAT 



iN-j){N 
is fulfilled. 

Proof: It is easy to prove the validity of the relation (|30D by simple straightforward 
calculations. □ 

Lemma 3: The matrix A{t) and vector d(r) as functions of r satisfy differential 
equations: 

A-\t)A'{t) = A'{0), A-\r)d'{T) = d'(r). (31) 

Proof: Taking into account group identities A{t + Ti) = A{r)A{Ti) and d(r + ri) = 
y4(r)d(ri) + d(r), we easy obtain 

A'{r) = hm + = A{r) lim ^ii^ = A{r)A'{0), 

Ti—>o n ri^o n 

d'(r) = lim + = A{r) lim ^ = A(r)d'(0) 

/^From these relations, it follows Eq. (|3TD . □ 
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Now, using differential invariants (^), we can construct invariant solutions. By turns, 
we define more convenient variables Uj by 



t/„ = w„, Uj = Vj-Y,{-lf(5j^kUj+kUt J=n-l,...,l. 

k=l 



Un = Un+nT = Un + C'^'^T. 



Observe that 
Next 

= Un-l + 2n/?„,„_irf/„ + n^/?„,„_ir^ 

Using the relations 
we obtain 

Thus, it is naturally to conjecture that the relation 



+ E Ci^lir'U^+k + Cr 1 (32) 

k=l 

holds for i = 1, ...,n. We observe that, consequently (H) and (|^), the relation (^) can be 
rewritten in vector form as 

u = A(r)U + d(r). 
Thus, we comes to the following ansatz for ^7v,n-invariant solutions 

u(x, t2, ...,tn) = A(r)U(X, T2, ...,Tn-2, Tn) + d(r). (33) 

In what follows, we will consider the case > 3. Taking into account Eq. ( pH] ) and lemma 
2, we obtain 



Proposition 2: Substitution of the ansatz (jj^ into the equations form = 1, A^— 
2, yields 

Ur^ = K„[U], m = l,...,N-2. (34) 
Here Km[U] are vector fields Km[u], where variables replaced by U-x ■ 



Proposition 3: Substitution of the ansatz ( pSj) into the system 
by virtue (p^, yields following constraint: 



'n \ ^ ^ fn 



i=l 



+ N-1] T^K;v-i[U] + E U + Ti+iM^] + A'(0)U + d'(0) = (35) 
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Proof: We have 



^ A'{r)V - Y ^ d'ir) 



[^ + N-l)tN 



(f + iV - l) 



+A{t) jux f-ai,i^ + 2a - ... + {-if-^N - 3)a,,j,J-^^^^ 



t 



N 

N~2\ 



J. j.N-5 



t 



N 



t 



N-3 
N , 



tN 



I = A{r) {k^^i[U] + 6^_i,jv-2rK^-2[U] + ... + 6^_i,ir^- 



Taking into account Eq. (|3^), lemma 3 and the relations 



n \ fn 

(iV - i) ( ^ + - 1^ ai,7v-i-i = + ^) «i+i,7v-i-2 + 



bN-i,i 



(f + iV-l 



N-i-2 



for A; = 2, - i - 2, z = 1, - 2, we comes to (|3|). □ 



Proposition 4: Substitution of the ansatz ^3[ ) into the system 
by virtue Eqs. and yields 

Ut, = K;v[U] 



Proof: We have 



{l + N-l) t 



(i + Ar-l)t^ 



-d'(T) 



+A{t) Vx ai, 



htN-i 

(-AT 



- 2ai,; 



(-AT 



+ ...-(-l)^-3(Ar-3)ai, 



tN-^t 



N-3 



N-2'^N-l 
N~3 



TV 



N-l 



1^2,2- 



f +2 

''4''Af-l 

t% 
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, ,Af-4 



(-l)--(iV - A)a,,^.:-^^^^ - {-ir-\N - 3)«2,iv-3^ ) + 



N-1 



N 



Now taking into account Eq. (p4D, lemma 3, the relations ( pq ) and the relation 

bN,i 



(f + iV-l 



we obtain (|37D. □ 



We observe that constraint (|35D is Galilean self-similarity condition 



- Q^,n(U) = E Q + T,+iK,[U] + A'(0)U + d'(0) = 0, 
where T]\f_i = 0. Consequence of (pH]) is scaling self-similarity condition 



i=l 



- Q^!n(U) = E ( ^ + ^ - 1 ) ^«K,[U] + JU = 0, 



where TAr„i = 0. By virtue (|^) and (|37D , we can replace in Eq. 

K,[U]^Ut„ 

for i = l,...,iV, i ^ N — 1, to obtain the constraint 



iV-ljT^UT^+ E (^- + ^-ljTiUT, + JU = 0. 



Example. Separately consider the case = 3. The ansatz 



u{x,h,h) = A{t)V{X,Ts) + d{T), T 



t2 



gives 



(f + 2)t3^ 



X = x- 



2(n + 4)U 



n 



- + 2 ) TgKaiU] + A'(0)U + d'(0) = 0, 



^ + 2j T,Ut, + ^XVx + JU = 0. 



More explicitly (^) is given by 



In (^^Uxxx + ^UUx 1 + 1 = 



15 



for n = 1 and 



n 



1. 



1 



n 



- + 2 ) T3 ( -UnXXX + UiUnX + T^UnUlX ] + U2 = 

+ 2 ) T3 ( -Ui-i,x + UiUnx + l:UnU,x ) + ^Ui+l = 0, z = 2, n - 1 



(43) 



-Un-l,X + -UnUnX ]+n = Q 



for n > 2. Eq. (|41|) , which presents scahng self-similar constraint can be immediately 
solved to yield 

f/.(X,T3) = ^/.(Z), Z = ^ 



(44) 



where 



P 



2{n-i + 1) 



n 



n + A n + 4 

Introducing (^4)) into (|4^) and (^31) we gives, respectively, 

^ {^Jzzz + ^//z) + 1 = 



(45) 



for n = 1 and 



1 



+ 2 -fnZZZ + /l/nZ + -/n/lZ + /2 = 



2+^ 



-/i-i,z + fifnz + -/n/iz ) + ifi+1 = 0, i = 2, n - 1 



(46) 



n 



-fn-l,Z + -fnfnZ ] + U = 0. 



for n > 2. 

According to Ablowitz — Ramani — Segur conjecture*, equations (^5]) and (|46|) should 
have the Painleve property. Single integration of Eq. ( ^5]) gives first Painleve trancedent^: 

f" + lf + lz = o 

whose canonical form is /" = 6/^ + Z. The system (|^) in the case n = 2 is reduced to 
second Painleve trancedent. Putting 

/l = ^/2 + ^Z, f2 = f 

and once integrating, we obtain second Painleve equation: 

f" + 2f + ^Zf + a = 

whose canonical form is /" = 2/^ + Zf + a. 

The problem to be addressed: to investigate Painleve property for the system (^), for 
every n. 



16 



^ L.V. Ovsiannikov, Group Analysis of Differential Equations, (Academic, New York, 1982). 

^ N.H. Ibragimov, Transformation Groups Applied to Mathematical Physics, (Reidel, Dor- 
drecht, 1985) 

^ P.J. Olver, Applications of Lie Groups to Differential Equations, Graduate texts in Math- 
ematics ( Springer- Verlag, New York, 1986). 

^ Y. Kosmann-Schwarzbach, Lett. Math. Phys. 38, 421 (1996). 

5 B. Fuchssteiner, Nonhnear Anal Theor. Meth. Appl. 3, 849 (1979). 

^ L.M. Alonso, J. Math. Phys. 21, 2342 (1980). 

^ A.K. Svinin, J. Phys. A: Math. Gen. 32, 5499 (1999). 

« M.J. Ablowitz, A. Ramani, H.J. Segur, J. Math. Phys. 21, 715 (1980). 

^ E.L Ince, Ordinary Differential Equations, (Dower, New York, 1956). 



17 



